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Abstract — In a recent work, a capacity-achieving scheme for 
the common-message two-user MIMO broadcast channel, based 
on single-stream coding and decoding, was described. This was 
obtained via a novel joint unitary triangularization which is 
applied to the corresponding channel matrices. In this work, 
the triangularization is generalized, to any (finite) number of 
matrices, allowing multi-user applications. To that end, multiple 
channel uses are jointly treated, in a manner reminiscent of space- 
time coding. As opposed to the two-user case, in the general case 
there does not always exist a perfect (capacity-achieving) solution. 
However, a nearly optimal scheme (with vanishing loss in the limit 
of large blocks) always exists. Common-message broadcasting 
is but one example of communication networks with MIMO 
links which can be solved using an approach coined "Network 
Modulation"; the extension beyond two links carries over to these 
problems. 

I. Introduction 

Multiple-input multiple-output (MIMO) Gaussian channels 
are a basic building block of many communication networks, 
due to their potential to enhance the throughput of communica- 
tion systems, and have been extensively studied both in terms 
of the theoretical limits (see, e.g., Hj) as well as in terms 
of modulation and coding schemes that allow to approach 
these limits. In different communication scenarios, different 
assumptions on the channel behavior and of the availability 
of channel state information are appropriate (see ||2], ||3] and 
references therein). 

A recent approach, coined "Network Modulation" @, 
tackles the problem of conveying information over different 
multiple-antenna multi-terminal networks where full channel 
state information is available at all terminals (i.e., a fully 
closed-loop scenario). The approach is based on jointly tri- 
angularizing several matrices using the same unitary matrix 
on one side (joint encoder or decoder) and different unitary 
matrices on the other side (separate decoders or encoders), 
such that the diagonals of the resulting triangular matrices 
satisfy desirable properties, e.g., that they are equal. This 
decomposition, along with successive interference cancellation 
(SIC) or dirty paper coding (DPC) ||5], transforms the channels 
into parallel scalar additive white Gaussian noise channels 
(AWGN). Thus, employing this scheme along with (any) 
scalar codes which are good for the AWGN channel, provides 
"practical" capacity-achieving schemes, for scenarios in which 
the capacity is known. Furthermore, somewhat surprisingly. 



it has been demonstrated that the approach allows to obtain 
new achievable rate regions to several information-theoretic 
problems, such as the two-way MIMO relay problem ||6l and 
the problem of joint source-channel coding of a source over a 
MIMO broadcast channel 101. 

A scenario of significant importance is that of sending a 
common message over a MIMO Gaussian broadcast (BC) 
channel, henceforth the multicasting scenario. The channel is 
given by 



1,2, 



(1) 



* This work was supported in part by the U.S. 
Foundation under grant 2008/455. 



Israel Binational Science 



where x is the complex-valued channel input vector of length 
n subject to a power constraint, (i = 1,2) is the output 
vector of user i of length m^. Hi is the rn^ x 7i complex channel 
matrix to user i, and is an additive circularly-symmetric 
complex Gaussian noise vector of length m^. Without loss 
of generaUty, we assume that both the noise elements and 
the input signal have unit power, i.e., Zi ~ 63^(0, Imt) 
and E [x^'x] < 1, where f denotes the conjugate transpose 
operator It was shown in [31 that it is possible to jointly 
triangularize the channel matrices Hi and H2 using unitary 
matrices, such that the ratio between the resulting diagonals 
is constant. This in turn allows to achieve the common- 
message capacity using single-stream encoding and decoding 
of standard AWGN codes along with SIC (much like in V- 
BLAST transmission for a single user ||7]). As we recall in 
the sequel, the problem of multicasting over several MIMO 
channels is tightly connected to the problems of universal 
coding over parallel channels, as well as rateless coding for 
Gaussian channels. Thus, the results we derive are relevant 
also to the latter problems. 

The joint triangularization of l|4j| was limited to only two 
matrices, and hence, only two-user multicasting (or "perfect 
two-rate" in the rateless problem [(8l) could be treated. The aim 
of the current work is to generalize the network modulation 
approach to more than two users. This is done by utilizing 
multiple uses of the channel, reminiscent of space-time coding 
techniques ll9l. ifTOl. 

II. Background: Network Modulation 

In this section we recall the joint triangularization of two 
matrices |4], and its application to the two-user multicasting 
problem. We then demonstrate the relevance of the scheme to 
the special case of a two-rate scalar Gaussian rateless problem. 



A. Unitary Matrix Triangularization 

The network modulation approach is based on several forms 
of matrix decompositions, one of which is the geometric 
mean decomposition (GMD) ITTl . For simplicity, we will only 
consider the decomposition of square matrices throughout this 
work. As we show in the sequel, this does not pose any 
restriction on the communication problems addressed. The 
GMD ifTTl of an ?7 x n matrix A is given by: 



A = UTV"^ , 



(2) 



where U and V aie nxn unitary matrices, and T is an upper- 
triangular n X n matrix such that all its diagonal values equal 
A, where A is a real-valued non-negative number 

Building on the GMD, the following decomposition, which 
will be referred to as Joint Equi-diagonal Triangularization 
(JET), was introduced in j4|. Let Ai and A2 be two complex 
matrices of dimensions n x n such that dct(A) = dct(_B). 
Then, the joint triangularization of Ai and A2 is given by: 



Ai = UiRiV'' 

A2 = U2R2V^ 



(3) 



where Ui,U2,V aie n x n unitary matrices, and i?i , R2 are 
upper-triangular nxn matrices with the same real-valued, 
non-negative diagonal values, namely. 



B. MIMO Multicast Scheme 



Vi = 1, 



We now recall how the JET decomposition can be used to 
obtain a practical scheme for transmitting a common message 
over a MIMO Gaussian BC with two receivers, as described by 
([T]i- Define the mutual information between a Gaussian input 
vector X, having a covariance matrix Cx = E [xx^] , and the 
channel output y^, by 



I{H,, Cx) = log dct (/ + H,C^Hj 



(4) 



The common-message capacity is given by the (worst-case) 
compound-channel capacity expression (see, e.g., lfT2l ): 



C 



max mill liHi, Cx) . 

Cx:tr(Cx)<l 4=1,2 



(5) 



Let Cx be an admissible covariance matrix, and assume for 
simplicity that I{Hi,C^) = /(i?2,Cx) = R. The following 
scheme IH achieves the rate R. 

Define the following augmented matrices: 



In 



where Fi = Hi^/CZ and /„ is the n x ji identity matrix. 

Next, the matrices Gi are transformed into square matrices, 
by means of the QR decomposition: 



Gi — QiGi 



(6) 



where Qi is an (nii+n) x n matrix with orthonormal columns 
and Gi is an n X n upper-triangular matrix with real-valued 



positive diagonal elements. Now, assuming that /(iJi,Cx) = 
I{H2,Gx) = R, this impUes ||4] Proposition 1]: 

det(Gi) = det(G2) = 2* . 
Therefore, Gi , G2 can be jointly triangularized using the JET: 

G, = U,R,V\ ^ = l,2, (7) 

where i?i and R2 are upper-triangular, having the same 
diagonal elements. The transmission scheme is as follows: 

1) Construct n optimal codes for scalar AWGN channels. 
The A;-th codebook is designed for a SISO AWGN chan- 
nel with a rate 21ogrfc, where is the k-th diagonal 
element of Ri (and also of i?2)- 

2) In each channel use, an n-length vector x is formed 
using one sample from each codebook. The transmitted 
vector X is then obtained using the following precoder: 



3) At the receiving ends, the z-th user calculates 



(8) 



(9) 



where Qi consists of the first n rows of Qi. 
4) Finally, the codebooks are decoded using SIC, starting 
from the n-th codeword and ending with the first one: 
The ji-th codeword is decoded first, using the n-th 
element of y^, treating the other codewords as AWGN. 
The effect of the n-th element of x is then subtracted out 
from the remaining elements of y. Next, the (n — l)-th 
codeword is decoded, using the (n — l)-th element of 
Yi - and so forth. 

The optimality of this scheme was proved in ||4] Sec. IV]. 
Example 1 (Application to the two-rate rateless problem): 

Consider the scalar Gaussian rateless problem defined in llTJI : 



Vn 



TO = 1,2, . 



The gain a is known only to the receiver, and can take one 
of M possible values, such that a gain of a,„ implies that the 
message should be decodable using m received blocks: Q 

i? = mlog(l + |a,„|^) , m = l,2,...,M. 

Specializing the problem to the case of one (possible) in- 
cremental redundancy block (AI — 2), the perfect two-rate 
rateless problem can be viewed as a 2-user MIMO-BC channel 
with channel matrices 



Hi 



ai 




H2 



012 
a2 

Applying the scheme of Section III-BI yields the following 
precoding matrix [Sj: 

2«/4 



V 



1 



1 



2R/2 + 1 1^ 2^/4 -1 
which coincides with the result in ifTsl Section III]. 

'Alternatively, this can be viewed as a scheme that works for every value 
of a, but designed to be optimal only for M specific values. 



Erez, Trott and Wornell |fT3l also treated the case of M = 
L = 3, and found a condition for which a "perfect" scheme 
exists. In the sequel we will shed light on this condition. 

III. Joint Triangularization of Many Matrices 

In this section we extend the network modulation technique 
to a any finite number of users, using a recursion principle. 
Specifically, given K matrices Gi , . . . , Gk, we wish to find K 
matrices with orthonormal columns Ui, . . . ,Uk, and another 
such matrix V, such that the matrices Ri = U^GiV are upper- 
triangular, having equal diagonals. We shall refer to this as 
A'-matrix JET, or simply A' -JET. 

The proof of the existence of a JET decomposition of two 
matrices Gi and G2 EJ is based upon applying the GMD (|2]i to 
the single matrix GiGj ^. Similarly, we show in the following 
lemma that {K + 1)-JET is equivalent to simultaneous GMD 
of K matrices, which will be referred to as A'-GMD. 

Lemma 1: Let Gi, . . . , G/f+i h& n x n complex valued 
matrices with equal determinants, and define the K matrices: 



where: 



K+l ■ 



(10) 



Then, there exist A' + 1 matrices with orthonormal columns 
f/i, . . . , Uk+1, of dimensions n x m, such that 



1, 



(11) 



where {Ti} are upper-triangular with all diagonal entries equal 
to 1, if and only if there exists an n x m matrix V with 
orthonormal columns, such that 



1,. 



where {Ri} are upper-triangular with equal diagonals. 

Proof: See a constructive proof in Appendix |A] ■ 
Remark 1: Constructing matrices with constant diagonals 
could be advantageous in practice, as this corresponds to equal 
gains of all the resulting sub-channels, and hence enables to 
use the same (single) codebook over all of them. 

We are thus left with the task of performing AT-GMD to 
AT matrices. In Section |IV] we state sufficient and necessary 
conditions for the existence of the above decomposition for 
the special case of two real-valued 2x2 matrices. We will 
then, in Section|Vl present a different approach, involving joint 
triangularization of block-diagonal matrices, which enables a 
nearly-optimal network-modulation scheme, even when exact 
triangularization is not possible. 

IV. Exact Triangularization with constant 

DIAGONALS OF TWO REAL-VALUED 2x2 MATRICES 

We now provide a necessary and sufficient condition for the 
existence of 2-GMD for real-valued 2x2 matrices. 

Theorem 1 (2-GMD for 2x2 real-valued matrices): Let 
Ai and A2 be real-valued 2x2 matrices with determinants 
equal to 1. Apply (any) JET decomposition to them: 



1,2. 



(12) 



ri Xi 
r2 



Then, there exist three complex-valued 2x2 unitary matrices 

C/GMD^ [;GMD^ yGMD ^^^^ that: 



1 * 
1 



if and only if the following inequality is satisfied: 

( Xi + X2\^ ^ , X1X2 

r2 7: < ?'2 + 



ri - r2 



(13) 



Without loss of generality, we can assume that the solution is 
of the form: 

aGMD _ / *1 32 



(14) 



Proof: The proof is straightforward, and is given in 
Appendix [C] ■ 

Remark 2: Although this theorem is valid only for real- 
valued matrices Ai and A2, the resulting unitary matrices 
Ui, U2, and V are, in general, complex-valued. In Section [V-AI 
we bring a restatement of the theorem, which involves only 
real-valued orthogonal transformations. 

Remark 3: This theorem can be applied to the three-rate 
rateless problem defined in Section III-BI This yields a con- 
dition for the existence of a perfect scheme, namely, R < 
6 log ( '^t/^ ) ~ 8.331, as in ifTSl . The details are given in 




Appendix |D 



V. Space-Time Triangularization 

As indicated by Theorem [T] joint triangularization with 
constant diagonal values is not always possible. However, even 
when the condition for joint triangularization does not hold, 
we can still perform nearly-optimal network modulation, by 
utilizing multiple uses of the same channel realization. The 
idea of mixing the same symbols between multiple channel 
uses has much in common with Space-Time Codes iU, ifTOl . 

A. Restatement of Theorem |7] 

In order to introduce the space-time like structure, we start 
by a restatement of Theorem [T| 

Recall the two-user common-message broadcast MIMO 
channel ([T|i with two transmit antennas (n ~ 2), and a general 
number of antennas m; at each receiver We now utilize 
transmission in two consecutive time instances (as in iQ)- This 
is equivalent to sending extended symbols over the following 
extended channel: 



i = 1,2. 



The extended vectors X, Y, Z are composed of two "physical" 
input, output, and noise vectors, respectively, and IK^ is the 
{2mi) X 4 extended channel matrix defined as (i = 1,2) 



^The JET decomposition is, in general, not unique. 



J{i = \Hi\ 



(15) 



where \A\ denotes the Kronecker product In ^ A, viz. a 
block-diagonal matrix with blocks of A on its diagonal: 



/AO 
A 



V 



\ 





The power constraint now becomes E [X-^X] < 2. 

Let Cx be a covariance matrix satisfying trace (Cx) < 1, 
and define the augmented matrices Gi as in (|6]l. Following 
Lemma [H we define the two 2x2 matrices: 



Ai — G\G^ ^ 



A2 — G2Gr, 



Also define the following 4x4 extended matrices {i 



1,2): 
(16) 



Since the matrices Ai and A2 are real-valued matrices, we 
can obtain 2-GMD of the matrices Ai and A2 under the same 
conditions as in Theorem [T] such that all the involved unitary 
transformations become real-valued. Following Lemma[T] this 
yields a 3-JET of the three matrices Si, S27 S3: 

where 3?^ are upper triangular with equals diagonals. 

In particular, the complex precoding matrix \/GMD gjygjj i^y 
(fT4l i implies the following (real) orthogonal space-time block 
code structure of V^^^ HI: 



GMD 



a 


— c 


b 


d 


\ 


b 


d 


—a 


c 




c 


a 


d 


-b 




d 


b 


c 


a 


) 



The same scheme as in Section III-BI can now be employed, 
such that the two channel uses are effectively transformed into 
four scalar AWGN channels, having the same capacities for 
all three users. Note that the matrix Qi of (|9]l is replaced with 



its extended version. 



Q^ 



®2 



B. Nearly-Optimal 2-GMD 

We now show how to utilize a space-time structure in order 
to obtain nearly-optimal joint triangularization of two matrices, 
such that the resulting triangular matrices have a constant 
diagonal. This method will later be generalized to any number 
of matrices, using Lemma [T] The resulting scheme becomes 
asymptotically optimal for large values of N, where N is the 
number of channel uses assembled together for the purpose of 
joint decomposition. Note that the proposed scheme is nearly 
optimal for any two complex-valued channels Hi (and not 
restricted to real-valued matrices, in contrast to the perfect 
construction of Theorem [T]i. 

Theorem 2 (Nearly-Optimal 2-GMD): Let Ai and A2 be 
two complex-valued nx n matrices, and define the following 
nN X nN extended matrices: 



Ai — \ Ai\ 



(S>N 



L2, 



(17) 



Then there exist three nNxn(^N— (n — 1)) matrices Hi, 1X2, V 
with orthonormal columns, such that: 



U]A,V = 



/ 1 * 
1 





V 



i = 1,2. 



1 / 



By using this decomposition together with Lemma [T] the 
same scheme as in Section ITl-BI can be employed, such that the 
N channel uses are effectively transformed into n{N — n + 1) 
scalar AWGN channels. The sum of the capacities of these 
channels tends to the capacity of the original channel for large 
values of N, where the only loss comes from edge effects 
(truncating the extreme n{n — 1) elements). 

The full proof of the theorem is given in Appendix |E] The 
main idea of the proof is demonstrated by the proof for the 
2x2 case, presented next. 

Proof of Theorem ^for n = 2: We start by jointly 
triangularizing the matrices Ai and A2: 



{ur^^A^v 



JET 



''1 




a-, 
r2 



(18) 



where rir2 = 1. We now apply the decomposition ( fTSl l to 
each block separately, using: 



U 



JET 



which yields the matrices 



.JET 



A JET I 



JET 



ri Xi 
^\r2 



■ I 


1 
" 


ri \ Xr 
' r2 





V 

Note that the sub-matrix 
A = 



I ^2 
I 









5 


























fl 







r2 



(19) 



does not depend on i, and therefore it can be decomposed 
using the GMD ©, A = {jCMDj^i-yGMD^t^ ^j^^^.^ 

upper-triangular with only Is on the diagonal. We use this 
decomposition to construct a second transformation - only 
this time it is not be applied on each block separately, but 
rather "mixes" pairs of consecutive blocks, using: 



GMD 



/ • 

|"[^GMD 

\ • 



\ 



®(N-\) 

■ 



GMD 



/ • • 

[yGMDJ 

\ •• 



• 

®(Af-l) 

• 



Applying this transformation to (fT^t yields the following 
{2N — 2) X {2N — 2) upper-triangular matrix: 



/I" 1 * I 
10 II 



u.U,v = 





V 



; 1 

I 



where U, ^ Uf ^uGmd ^ a ^jet^gmd^ , 

C. Nearly-Optimal K-GMD 

By using Lemma [T] we can generalize Theorem |2] to any 
number of users, as follows: 

Theorem 3 (Nearly-Optimal K-GMD): Let Ai,. . . , Ak be 
K complex-valued nxn matrices with determinants equal to 
1, and define Ai,. . . ,Ak as in ( fTTj i. Then there exist K + 1 
matrices Hi, ... , VLk, V, with orthonormal columns, such that: 









V 



1^1,. ..,K. 



1 / 



Proof: A sketch of the proof is given in Appendix |B] ■ 
VI. Discussion 

Theorem [U provides sufficient and necessary conditions for 
joint GMD of two real-valued 2x2 matrices. This naturally 
raises the question of how this condition can be carried over 
to the complex-valued case, and to general dimensions nxn. 

Furthermore, we demonstrated that (exact) K-GMD, not 
using any space-time structure, is not always possible. Nev- 
ertheless, near/y-optimal communication schemes can always 
be constructed, which become optimal in the limit of large N. 
It remains an open question whether an exact triangularization 
can be obtained using only a finite number of channel uses. 
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Appendix A 
Proof of Lemma[T] 

Proof of LemmalJ} The direct part holds trivially. We are 
therefore left with the task of proving the converse part. We 
start with the QR decomposition G^^jC/^+i = VS, where 
V is of dimensions nxm with orthonormal columns, and S is 
an m X m upper-triangular matrix. Thus, using (ITOi and (ITll . 
we obtain the following equalities: 



UjG^VS = T, , i = 1, . 
uI^,Gk+iVS = I . 
Multiplying by S^^ on the right yields: 



,K 



u}G,V = T,S-^ , 
uI^^Gk+iV = 



1 = 1,. ..,K 



Since Ti are upper-triangular with only Is on the diagonal, the 
matrices = T,S~^ (i = 1, . . . , K) and Rk+i = 5"^ have 
equal diagonals, which completes the proof. ■ 

Appendix B 
Sketch of Proof of Theorem[3] 

Proof Idea: The theorem has already been proved for the 
special case of K = 2. For larger values of K we prove by 
induction, applying repeatedly Lemma [T| and of Theorem |2l 

1) According to Lemma [1] performing /\-GMD is equiv- 
alent to {K + 1)-JET. We can thus transform K upper- 
triangular matrices with constant diagonal values into 
K + 1 upper-triangular matrices of the same size, 
Ri, . . . , Rk+1 with equal diagonals. 

2) Given the matrices . . . , Rk+i, construct the block- 
diagonal extended matrices "Ri, as in ( [TtI i. Using the 
technique of Theorem |2] we construct matrices with or- 
thonormal columns, U^^^^^\. . . , U^'^++^\ V^-^+i), such 

that the matrices {11^^^^^%^^^+'^'^ are upper- 
triangular, with constant diagonals. Finally, the loss in 
rate could be made arbitrarily small by taking N to be 
sufficiently large. 

Appendix C 
Condition for 2-GMD of Real-Valued 2x2 
Matrices 

We now prove the necessary and sufficient condition for the 
existence of joint-triangularization of two 2x2 real-valued 
matrices. 

Proof of Theorem Q} Let Ai and A2 be real-valued 
2x2 matrices with determinants equal to 1. Apply the JET 
decomposition to these matrices, to obtain 



A = C/JETp.fyJETNt 



1,2, 



(20) 



where: 



Ri = 



ri Xi 
ra 



such that rir2 = 1. The matrices U^^^ , V'^^^ are real-valued 
unitary matrices, and we can assume without loss of generaUty 
that dot (yJET^ ^ ^_ 

If there exist three complex-valued 2x2 unitary matrices 
Ui,U2,V such that: 



(21) 



then according to (|20] |. 



GMD\t 



GMD 



1 * 
1 



where 



rrGMD 
^ i 

yGMD 



JET\ t 



V . 



Denote the entries of the first column of v'^^^ by si and Sj, 
i.e., 

yGMD ^ ^1 * 
V ^2 * 

where |sip + |s2p = 1. The first column of RiV'^^^ and of 
i?2T/G!MD therefore: 



R^V 



GMD 



1,2, 



riSi + XiS2 * 
'"2S2 * 

where xi and X2 denote the off-diagonal elements of Ri and 
i?2 respectively. These two columns must have a norm of 1, 
namely; 

|risi + XjSsI^ + k2S2|^ = 1 , i = l,2. 

Since ri,r2,x are real-valued, si and S2 must satisfy the 
following three equations: 

kiP + |.2p = 1 

rl\si\^ + {xl + rl)\s2\^ + 2riXiRe{siS2) = 1 
r?|sip + (x^+r|)|s2p + 2ria;2Re(siS2) = 1- 
Denoting a = and substituting a in these equations, results 



l+\a\' = 

rj|ap + (a;? +r^) + 2riXiRe(a) = 

rj|ap + (a;^ +r^) + 2riX2Re(a) = 

Subtracting (l24l i from ( |23] l yields: 

(xi +a;2) + 2riRc(a) = 0, 

So we have: 



1 



1 

N 
1 

N 



(22) 
(23) 
(24) 



S2 



- 1 



Rc(a) 



(Re(a))^ = 



(Im(a))^ 



Xi + X2 



Xi + X2 



S2 



r2 



Xi + X2 



(25) 
(26) 
(27) 
(28) 



Thus, equation (|23] | and ( l24l i become: 

- 1 



X1X2 



and therefore equation (|27i becomes 



(Im(«))^ 



2:1X2 



1 



1 



Xi + X2 



Thus, the following conditions are necessary and sufficient for 
the existence of a solution: 

rf - 1 



■ X1X2 



r| + 0:10:2 
r? - 1 



0:1 + 2:2 



> 



> 0, 



which are equivalent to 

'0:1 + 0:2 

f 2 ■ 



< r2 



X1X2 



(29) 



n - r2 

This proves that i29[ is a necessary condition for the 
existence of the decomposition ( I2TI ). 

Now, assume that this condition holds, and define the 
matrix: 

yGMB _ f Si S2 

We now apply the QR decomposition to the matrices RiV'''^^ 
and i?2'^"™°: 



/^GMDNtp.^GMD 







(30) 



The first columns of both RiV'^^^ and R2V'^^^ have 
norms equal to 1. Therefore, from the construction of the QR 
decomposition, it follows that ai = 02 = 1. Consequently, 
since both matrices have a unit determinant, ci = C2 = 1 
must hold as well. Thus, (l30b becomes: 



and therefore. 



GMD\ t 



)'R^V 



GMD 



A,. 



1 * 
1 



1 * 

1 



1/t, 



where 

= 

V = 

Furthermore, since the matrix I/GIMD j-j^g form 



rrJETrrGMD 
^ i ^ i 
^JET^GMD 



Si 

Sn 



32 

-sT 



yGMD 

^2 ■^l 

and V'^^'^ is a real-valued unitary matrix with unit determi- 
nant, it is easy to see that the matrix V is also of the form 



V 



51 S2 

52 —s-y 

which completes the proof of the theorem. 



Appendix D 
Three-Rate Rateless 

We now consider the three-rate "rateless" problem, as 
defined in Section III-BI with M = L = i and a given rate R: 



15^1 



IH2 = 



ai 













•) 










a2 





^ 





Ct2 











^ 


as 











as 











"3 



where ai,a2,a3 are the positive values satisfying log(l 



at 



21og(l + ^ 3iog(i 



§) = R. As in the 2- 
rate case, the covariance matrix in this problem is the identity 
matrix, Cx = /. Since J-C^ is a s scaled identity matrix, we 
can ignore it and concentrate on the remaining two matrices. 
The augmented matrices, as defined in (|6]l, are: 



Si = 



2f 











1 











1 


2- 










2f 













The decomposition (fTZt becomes; 



3?i 



3^2 



where 



2- 



1 




1 



2f. 



It then follows from Theorem [T] that there exists a perfect 
solution over the complex field if and only if: 

™2 



or explicitly: 

2- 



4 < 0, 



1 - 3 • 2- + 2- < 



1 + 2- 

This condition is satisfied if and only if: 

'3 + V^' 



i? < 6 log 



8.331 



which coincides with the result that was obtained in fT3]. 

For rates higher than this threshold, a perfect capacity- 
achieving solution does not exist. However, as explained 
earlier, multiple channel usages can be utilized in order to 
approach capacity asymptotically. 

Appendix E 
Nearly Optimal 2-GMD for n>2 

We now bring the proof of Theorem |2] for the general case 
n>2. 

Proof of Theorem |2} 

Let Ai and A2 be the two nx n complex valued matrices, 
with determinants equal to 1. As in (fTSl l. we define the 
extended matrices, 

(A, • • • \ 
A, ■■■ 



A, = 



\ 



We are looking for three nN xnlN ~n+l) matrices Hi, 1X2, V 
with orthonormal columns, such that 

n\A,v 

are upper-triangular, with only Is on the diagonal. 
We accomplish that using three steps: 

a) Joint Triangularization: As in the n = 2 proof, we 
start by jointly triangularizing the matrices Ai and A2: 



r2 



V • • • r„ / 
We now apply this transformation to each block separately: 



A: 



JET 



JET\ t J ^oJET 



(31) 



where 



U: 



JET 











JET 



v 









.. [/JET 
























. . l^JET 



\ 



b) Reordering: It will now be convenient to re-order the 
columns of A^^"^ such that the following columns: 

fcn, kn + [n — 1), kn + 2{n — 1), • • • , fcn + (n — 1)^ 

will become "grouped together" for every fcH Formally, We do 
so by introducing the following nNxn{N—n+l) re-ordering 
matrix 0: 

= { Otherwise, ^^^^ 

where the function tt is defined as follows: 
. For 1 < j < n{N -n + 1), 



TTj = (n - 1) • [{j - l)%n] + : 



.7-1 



As a result of this re-ordering, we obtain an upper- triangular 

[N — n + l)n X {N — + l)n matrix, which has a block- 
triangular structure: 



/ 4i * 1 

I A I 
- -I- - -I- 
: I : I 

1 4- 1 4- 
I I 



* I * \ 

* I ^ 

■ . I : 

- - 4- - 

A I * 



VOIOI---I 'A/ 
^11 I I ' 



Note that this set includes exactly one symbol from each of n consecutive 
channel uses. 



where 



A = 





Tn-l 





\ 





V ••• ri / 

c) GMD: Since the matrix A does not depend on i, we 
can decompose it using GMD: 

GMDji/'^GMD^t 



where T is upper-triangular with only Is on its diagonal. 
We now apply the GMD to all the blocks of O'^Aj^^O: 



w 



GMD 



V' 



GMD 



/ jjGMD 




V 

/ yGMD 







jjGMD 




yGMD 






jjGMD 






V 



GMD 



[6] A. Khina, Y. Kochman, and U. Erez, "Physical-layer MIMO relaying," 

in Proceedings of ISIT-201 1, St. Petersburg, Russia, July 2011. 
[7] P. W. Wolniansky, G. J. Foschini, G. D. Golden, and R. A. Valenzuela, 

"V-BLAST: an architecture for realizing very high data rates over the 

rich-scattering wireless channel," in ISSSE 1998, URSI International 

Symposium, pp. 295-300. 
[8] A. Khina, Y. Kochman, U. Erez, and G. W. Womell, "Incremental coding 

over MIMO channels," submitted to Infomiation Theory Workshop 

(1TW2011), Parati, Brazil. 
[9] S. Alamouti, "A simple transmit diversity technique for wireless commu- 
nications," Selected Areas in Communications, IEEE Journal on, vol. 16, 

no. 8, pp. 1451 -1458, Oct. 1998. 
[10] V. Tarokh, N. Sheshadri, and A. R. Calderbank, "Space-time block codes 

form orthogonal designs," IEEE Trans. Information Theory, vol. IT-45, 

pp. 1456-1467, July 1999. 
[11] Y. Jiang, W. Hager, and J. Li, "The geometric mean decompostion," Lin. 

Algebra and Its Applications, vol. vol. 396, pp. 373-384, Feb. 2005. 
[12] D. Blackwell, L. Breiman, and A. J. Thomasian, "The capacity of a 

class of channels," The Annals of Mathematical Statistics, vol. 30, pp. 

1229-1241, Dec. 1959. 
[13] U. Erez, M. D. Trott, and G. W. Womell, "Rateless coding for Gaussian 

channels," 2007. [Online]. Available: http://arxiv.org/abs/0708.2575 
[14] A. V. Geramita and J. Sebeny, "Orthogonal designs, quadratic forms and 

Hadamard matrices," Lec. Notes in Pure and App. Mathematics, 1979. 



to obtain: 



^'^GMD'^tQTgJETQ'^GMD 



I 



I 7"; I • • • I * I * 
1 I 1 I 

: I : I : i ■ . i : 
I- i-i-lui 

V I I • • • I I 7 

' ' ' ' (33) 
where Ti are upper-triangular with only Is on the diagonal. 
We conclude by combining ( l33b with dSTI ) to obtain: 




V I 



where 



V 



yJETQ^GMD 
-^JETq-^GMD 



which completes the proof of the theorem. 
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